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One of the pivotal questions in the physics of high-temperature superconductors is 
whether the low-energy dynamics of the charge carriers is mediated by bosons with a 
characteristic timescale. This issue has remained elusive since electronic correlations 
are expected to dramatically speed up the electron-boson scattering processes, confin¬ 
ing them to the very femtosecond timescale that is hard to access even with state-of- 
the-art ultrafast techniques. Here we simultaneously push the time resolution and the 
frequency range of transient reflectivity measurements up to an unprecedented level 
that enables us to directly observe the ~16 fs build-up of the effective electron-boson 
interaction in hole-doped copper oxides. This extremely fast timescale is in agreement 
with numerical calculations based on the t-J model and the repulsive Hubbard model, 
in which the relaxation of the photo-excited charges is achieved via inelastic scattering 
with short-range antiferromagnetic excitations. 


After almost 30 years of intensive experimental and theoretical efforts to understand the origin of high-temperature 
superconductivity in copper oxides, a consensus about the microscopic process responsible for the superconducting 
pairing is still lacking. The large Goulomb repulsion U eV between two electrons occupying the same lattice 
site is believed to have fundamental consequences for the normal state of these systems [T], and it is not clear 
whether a BGS-like bosonic glue that mediates the electron interactions and eventually leads to pairing can still be 
defined [ng. The fundamental issue can be reduced to the question whether the electronic interactions are essentially 
instantaneous, or whether the low-energy physics, including superconductivity, can be effectively described in terms of 
interactions among the fermionic charge carriers mediated by the exchange of bosons. The problem can be rationalized 
by considering the Hubbard model, in which the instantaneous virtual hopping of holes into already occupied sites 
(with an energy cost of U) inherently favors an antiferromagnetic (AF) coupling J=At\/U between neighbouring 
sites, where tf, is the nearest-neighbour hopping energy. As a consequence, antiferromagnetic fluctuations with a 
high-energy cutoff of 2J<^U naturally emerge as a candidate for mediating the low-energy electronic interactions, on 
a characteristic retarded timescale on the order of h/2J. 

In principle, time-resolved optical spectroscopy [5] may be used to prove the existence of an effective retarded boson- 
mediated interaction, provided that: i) the temporal resolution is on the order of the inverse bosonic fluctuation scale 
(e.g. h/2J for AF fluctuations); ii) the optical properties are probed over a sufficiently broad frequency range, to 
obtain direct information on the dynamics of the electron-boson coupling. Here we develop a transient reflectivity 
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FIG. 1: Equilibrium and non-equilibrium optical properties of YBi2212. a) The normal state reflectivity (black line) of the 
YBi2212 crystal is reported. The dashed line represents the reflectivity calculated by increasing the electron-boson scattering 
rate, 700 , in the extended Drude model, b) The ultrafast dynamics of 5R/R{u),t), measured across the isosbestic point u), is 
reported. The measurements have been performed with an incident pump fluence of 0.7 mj/cm^. The linearity of the optical 
response has been checked in the 0.2-2 mj/cm^ fluence range, c) The first eigenfunction, obtained through the SVD (see 

Supplementary Information) is displayed (yellow dots). SVD consists of decomposing 5J?(a;,t)/i?eq(ti;)~5])j,'kfe(a;)(5(()fc(t). The 
fit obtained by increasing the total scattering rate, i.e., i(u)) = l/Req(oj)-dR /is reported as a black line, d) The time- 

dependent scattering rate variation, 57 oo(t)= 5 <()i(t), is reported (blue dots). The black line is the fit to the data of the function 
(l-exp(-t/rr))-(Iiexp(-t/ri)-|-l 2 exp(-t/r 2 )) convoluted with a gaussian pulse accounting for the experimental time-resolution. 
The red line is the measured cross-correlation trace between the pump and probe pulses, that sets the time-resolution to 19±2 
fs. The red arrow highlights the time delay (~40 fs) at which the maximal 57cx) variation is measured. 


experiment |S], based on extremely short broadband pulses (9-13 fs) in the infrared-visible spectral region (0.75-2.4 eV) 
(see Methods), that is most sensitive to the variations of the electron-boson scattering rate. The goal is to investigate 
possible retardation effects of the electron-boson interaction that were inaccessible in previous measurements with 
~100 fs time resolution [7]. The qualitative idea of the experiment is to use an ultrashort light pulse to impulsively 
increase the kinetic energy of a small fraction of the charge carriers (photo-excited holes) in doped cuprates. A second 
broadband pulse probes the instantaneous optical scattering rate by inducing boson-assisted optical transitions of 
the charge carriers in the conduction band. As long as the number of photo-excited holes is small as compared to 
the intrinsic doping of the system (see Methods), we can assume that we are probing the average scattering rate 
experienced by the holes which have not been directly excited by the pump pulse. Since the optical scattering rate is 
proportional to the energy stored in the bosonic fluctuations, we can monitor the time needed to complete the energy 
exchange from the photo-excited carriers to the bosonic bath. 

The dynamics of the optical scattering rate after the impulsive excitation can be directly monitored by probing the 
damping of the infrared plasma edge of a prototypical superconducting cuprate. This is easily shown by inspection 
of the room-temperature equilibrium reflectivity, Req{uj), of a Bi2Sr2Yo,o8Cao.92Cu208+5 (YBi22I2) crystal [8]. The 
doping level (hole concentration p~0.13, Tc~83 K) is close to optimal (p~0.16, Tc~96 K), and we estimate the 
opening of the pseudogap to occur at T*~240 K [ 3 . As shown in Figure la, the infrared region is dominated by 
a broad plasma edge, which can be effectively reproduced by an Extended Drude Model (EDM) [10] described by 
47rcr£i(a;)=a;p/[7(a;,r)-ia;(l-|-A(a;,r))] where aoi^j) is the EDM optical conductivity, uip the plasma frequency and 
1-l-A the renormalization factor of the effective mass of the carriers. The optical scattering rate, 7 ( 0 ;, T), controls the 
damping of the plasma edge and it is, in principle, a function of the frequency and temperature, as a consequence 
of the coupling to bosonic fluctuations whose occupation number is regulated by the Bose-Einstein statistics m- 
Nevertheless, in the energy window probed in the experiment, 7 ( 0 ;, T) varies less than 10% as a function of frequency 
(see Supplementary Information), therefore we can safely discuss the time-resolved results considering its asymptotic 
(w—too) value, i.e., 7 oo(T). If the total scattering rate is enhanced (( 57 oo> 0 ) by any effect, this leads to a further 
damping of the edge and to the evolution of Req{uj) shown in Eigure la. An isosbestic point is found at the frequency 
a;~l.l eV, across which the reflectivity variation changes from negative to positive. 
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FIG. 2: Microscopic calculation of the energy transfer from a photo-excited hole to local antiferromagnetic bonds through 
the out-of-equilibrium t-J model (see Methods for further details), a) The time evolution of the hole kinetic energy variation, 
Ai7kin(t). In the main panel we set J/th=0-S, where th=360 meV represents the hopping amplitude. The inset displays 
the range of the characteristic relaxation time Ttj at different values of J/th, after a finite-size scaling (see Supplementary 
Information), b) Snapshots of the relative energy increase of the local antiferromagnetic bonds Ai?a(t,r) on a 9x9 square 
lattice around the photo-excited hole (red circle). The thickness of the black segments is proportional to the energy stored in 
each bond. The grey dots in a) indicate delay times of the snapshots. 


In Figure Ib we report the dynamics of the relative reflectivity variation, 6R{u!,t)/Req{uj), measured at T=300 K 
on YBi22I2 as a function of the delay t between the pump and probe pulses. The sub-10 fs pump pulse is set to 
2.1 eV central energy, in order to avoid the spurious signals measured when pumping below the Drude absorption 
edge m- The probe spans the energy range across oj, that is more sensitive to the possible variations of 'yoo{T), 
as already shown in Fig. la. One key result of this experiment is that, for each given time delay t from the 
very first femtoseconds up to the picosecond timescale, 5R(uj,i)/Req(oj) can be accurately reproduced exclusively by 
increasing 700 . This is corroborated by the quantitative analysis through Singular Value Decomposition (SVD) (see 
Supplementary Information), confirming that 97% of the signal is reproduced by factorizing the reflectivity variation 
in 5R{u),i)/Req{oj)=l/Req{uj)-dR /where ( 57 oo(t) contains the dynamics of the scattering rate variation. 

Remarkably, the maximum variation of 1^700 (t) is measured at a finite delay time (~40 fs) from the excitation, as 
shown in Figure Id. The intrinsic build-up time of the d 7 oo(t) signal (Tr=16±3 fs) is extracted by fitting the dynamics 
with the simple function (l-exp(-t/Tr))-(Iiexp(-t/ri)+l 2 exp(-t/r 2 )) convoluted with a gaussian pulse accounting for 
the experimental resolution (pump-probe cross-correlation with a full width at half maximum of 19±2 fs), that has been 
carefully checked through the Cross Frequency-Resolved Optical Gating method (see Supplementary Information). 
This result demonstrates that the energy exchange between the photo-excited holes and the boson bath takes a non¬ 
zero time of the order of 15-20 femtoseconds. The two further decay scales (ti~90 fs and r 2 "^l ps), measured in the 
6R{ui,t) signal, are related to the subsequent coupling with the optical buckling and breathing phonons (ti) and with 
the rest of the lattice vibrations {T 2 ) [7l IT^ IT^ . 

The present high time-resolution transient reflectivity experiment demonstrates that indeed the electron dynamics 
in the doped cuprates can be captured by an effective model in which the charge carriers exchange energy with boson 
excitations on an very fast (but non-zero) timescale. Considering that the intrinsic time fs) associated with 

the Cu-O-Cu hole hopping process inevitably leads to the creation of local spin excitations, AF fluctuations are a 
candidate as the dominant mediators of the hole interactions on this ultrafast timescale. As a further support of our 
assignment, neutron and x-ray scattering experiments [IlHIS] have revealed a rich magnetic dynamics in hole-doped 
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FIG. 3: Ultrafast dynamics of the electron-boson interaction. The different relaxation dynamics in high-Tc superconductor and 
in conventional superconductors are sketched in the left panel, a) The dynamics of (t) is shown for different superconducting 
systems (YBi2212, blue circles; BSLCO, red circles; od-BSLCO, green circles; MgB 2 , grey squares). The red arrows highlight 
the time delay corresponding to the maximal value of Sy^o (t) after the excitation, b) The critical temperatures of the systems 
considered is reported as a function of the build-up time (tt) of 5yoo{t), measured through the time-resolved optical spectroscopy. 
The black lines are the critical temperatures calculated by the Eliashberg-like expression Tc=0.83flexp[-1.04(l-|-A6)]/[Ai,-/r*(l -I- 
0.62A6)] |26| . fixing /r*=0 and assuming a generic relation between the electron-boson energy exchange rate, l/r^, in cuprates 
and the total coupling constant A;, and 17=200 meV (solid line). The dashed line represents the critical temperatures calculated 
for 17=70 meV and Ai,=l.l compatible with the phonon-mediated pairing in MgB 2 m- 


cuprates, characterized by resonance modes at 50-60 meV and a very broad spectrum with a cutoff of the order of the 
bandwidth 2 J~240 meV found for the AF parent compounds. Similar bosonic-fluctuation spectra have been extracted 
from equilibrium optical spectroscopy [lO]. On the other hand, we can rule out a possible photo-induced coherent 
rearrangement of the oxygen atoms, that has been demonstrated on the 20 fs timescale m, since: i) no significant 
variation of the charge transfer excitations (~2 eV), that are the most sensitive to a rearrangement of the oxygen 
orbitals, is observed (see Supplementary Information); ii) no coherent oscillation is detected in the 6R{u!,t)/Req{uj) 
signal at any wavelength. Even though a single incoherent scattering process with optical phonons could be extremely 
fast (1/4 of the period, i.e., ^20 fs for the highest energy modes at 70 meV), several scattering processes are necessary 
to complete the energy transfer to the phonons, as will be shown later in the case of the conventional superconductor 
MgB2. 

The role of short-range AF correlations as the fastest effective bosons mediating the charge dynamics is corroborated 
by studying the non-equilibrium dynamics of the photo-excited holes in the t-J Hamiltonian MM (see Methods), 
which is the minimal model that fully retains the dynamics at the energy scale J. Since the rate of the photo¬ 
excitation of holes in the system is very low (see Methods), we conjecture that the key relaxation mechanism at very 
short times corresponds to the creation of high-energy AF excitations in the close proximity of the photo-excited 
holes. We model the photo-excitation of a hole immersed in the AF background by instantaneously raising its kinetic 
energy, i.e., Ai5kin(t)=2 eV. The relaxation dynamics is calculated by numerically integrating the time-dependent 
Schrodinger equation, without any assumption of quasi-thermal intermediate states, as usually described by effective 
electronic and bosonic temperatures. The dynamics of AAkin(t) after the quench is reported in Fig. 2a. While the 
total energy remains constant during the time evolution, the observed relaxation of A Akin (t) is compensated by the 
simultaneous increase of the energy stored in the AF background. Figure 2b displays snapshots of the excess spin 
energy density, AAs(t,r), in the vicinity of the photo-excited hole. We note that the excitation process involves only 
the neighbouring bonds, demonstrating that short-range AF correlations are sufficient for describing the relaxation 
dynamics. The characteristic relaxation time, r^j, of the energy transfer to AF fluctuations is governed by the ratio 
J/th, as shown in the inset of Fig. 2a. Setting J/th=0.3 to a realistic value and performing a finite-size scaling 
analysis (see Supplementary Information), we obtain fs, that is in notable quantitative agreement with the 

experimental observations. We remark that, even though the t-J model cannot account for the optical properties at 
the energy scale U, the description of the low-energy electron dynamics through an effective charge-boson coupling 
is also supported by the generic single-band Hubbard model at equilibrium. In the Supplementary Material we use 
Dynamical Mean Field Theory (DMFT) to calculate the temperature-dependent optical conductivity of a hole-doped 
system (p=0.16). We demonstrate that the effect of an energy increase (in the form of a 5T variation) can be 
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rationalized as an increase of the number of bosonic fluctuations at the energy scale J without any significant change 
of spectral weight at the scale U. 

Finally, valuable insight into the general relation between the electron-boson coupling and the dynamics of <57oo(t) is 
provided by the results obtained on different families of cuprates and on more conventional superconductors (see Fig. 
3). We repeated the same experiment on optimally doped (BSLCO, a;=0.4, Tc=33 K) and overdoped (od-BSLCO, 
x=0.2, Tc=19 K) single Cu-0 layer Bi2Sr2_xLaa;Cu06 crystals. The time necessary for attaining the maximal variation 
of 7oo is non-zero (rr=21±4 fs for BSLCO and rr=22±4 fs for od-BSLCO), slightly longer than for YBi2212. These 
results demonstrate that the coupling to bosonic fluctuations on the ultrafast timescale is a general property of the 
normal state of hole-doped cuprates even at large hole-doping concentrations. Qualitatively different results are 
obtained on MgB 2 , that is considered the conventional phonon-mediated superconducting system with the highest 
critical temperature (Tc=39 K). Here, the maximum variation of the scattering rate after the impulsive excitation 
is significantly delayed (^90 fs, as shown in Fig. 3a) and rj.=32±3 fs is measured. This result demonstrates that, 
in the systems in which the superconducting pairing is mediated by high-energy meV) optical phonons [11) . 

the electron-phonon energy exchange process is completed on a timescale that is few times the inverse energy scale of 
the phonons involved. Furthermore, a single slower exponential decay (ti=500 fs) is measured in the 6R{uj,t) signal 
of MgB 2 , corresponding to the energy exchange with the rest of the lattice. On conventional non-superconducting 
metals, like Cu and Au, time-resolved measurements have been widely applied [3 [25], evidencing an electron-phonon 
coupling on the order of ^600 fs in Au and ^1 ps in Cu. Fig. 3 summarizes all the results obtained in both 
conventional and unconventional superconductors, demonstrating a universal relation between the energy scale of 
the boson modes mediating the electron interactions in the normal state and the timescale of the dynamics of the 
electron-boson scattering rate after the pump excitation. Intriguingly, the dependence of the critical temperature 
of the system on the experimental r^. follows the general and qualitative trend of the approximate solution of the 
McMillan’s equation [^, once we assume a direct relation m between the measured 1/r,. and the total coupling 
constant Xb=2 Jn(r2)/Hdr2, where n(H) is the coupling function to antiferromagnetic fluctuations (/^x(f2)) for copper 
oxides and to phonons (a^F(H)) for MgB 2 and metals. 

Although we do not exclude that the unmediated charge-charge interactions could influence the relaxation of very 
high-energy excitations (i.e. those involving the doublon formation as a consequence of a charge-transfer transition), 
our results are consistent with a picture in which, close to the optimal hole concentration, the dynamics of the 
charge carriers in the conduction band can be effectively described in terms of the interaction with short-range 
antiferromagnetic fluctuations on a timescale (~16 fs) on the order of few times /i/2J. This extremely fast dynamics 
has been confirmed by very recent non-equilibrium DMFT calculations [28] . We stress that, even in this picture, U is 
the crucial interaction since it determines the dynamics through the strength of the spin-charge coupling, J=4d.f^/U. 

Our results represent a new time-domain benchmark for realistic microscopic models of the pairing in high- 
temperature superconductors and pave the way to unveiling the pairing mechanisms in a wealth of intriguing systems, 
like different families of cuprates, iron-based superconductors [23 [53 all the materials that are on the verge 
of a (ferro-)antiferro-magnetic phase transition |4|. Furthermore, we provide a platform for the next generation of 
time-resolved experiments aimed at addressing the role of charge-fluctuations, the dual itinerant and local character 
of the charge carriers that would require multi-band models, and the additional coupling to bosonic modes related to 
the pseudogap phenomenon. 
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Methods 


A. Experimental setup 

A Ti:sapphire amplifier (Clark-MXR model CPA-1) delivers a train of pulses at 1 kHz repetition rate with 150-fs 
duration at 780 nm central wavelength and is used to simultaneously drive three Non-collinear Optical Parametric 
Amplifiers (NOPAs) operating in different frequency intervals. All NOPAs are seeded by white light continuum (WLC) 
generated in a sapphire plate. The first NOPA (NOPAl), is pumped by the fundamental wavelength and is amplified 
by fulfilling quasi-phase matching condition in a periodically poled stoichiometric LiTaOa (PPSLT) crystal. The signal 
generated in this process covers a spectral range between 1 prm (1.24 eV) and 1.5 /rm (0.83 eV) and is temporally 
compressed to nearly the transform-limited (TL) pulse duration (8.5 fs) by a deformable mirror based pulse shaper 
m- The second NOPA (NOPA2) is pumped by the second harmonic and amplifies in a beta-barium borate (BBO) 
crystal pulses with a spectral content between 820 nm (1.5 eV) and 1050 nm (1.2 eV), which are compressed to nearly 
TL 13-fs duration by a couple of fused silica prisms. Both these NOPAs serve to probe the transient response of the 
sample and are synchronized with a third NOPA (NOPA3), pumped by the second harmonic and using BBO, which 
initiates the dynamics. The spectrum of NOPA3 spans a frequency range between 510 nm (2.4 eV) and 700 nm (1.8 
eV) and it is compressed to 7 fs duration by multiple bounces on a pair of chirped mirrors, making the overall temporal 
resolution of the pump-probe setup below 19 fs. The time delay between pump and probe is adjusted by a motorized 
delay stage and both the beams are focused on the sample by a spherical mirror in a quasi-collinear geometry. The 
spectrum of the NOPAl and NOPA2 probes is detected respectively by an InGaAs and a Si spectrometer working at 
the full 1 kHz laser repetition rate. By recording the reflected probe spectrum at different temporal delays t with and 
without pump excitation, we measure the differential reflectivity: SR{u!,t)/Req{uj)=[R{ui,t)-Req{uj)]/Req{ui). 


B. Samples 


The Y-substituted Bi2212 single crystals were grown [8] in an image furnace by the travelling solvent floating-zone 
technique with a non-zero Y content to maximize T^. The underdoped samples were annealed at 550 °C for 12 days in 
a vacuum-sealed glass ampoule with copper metal inside. La-Bi2201 crystals were all by the floating-zone technique, 
and characterized as described elsewhere [35]. T12201od single crystals were grown [33] by a copper-rich self-flux 
method, with stoichiometry Tli 88(i)Ba2Cui.ii(2)06-i-5 corresponding to Cu substitution on the T1 site, thus away 
from the Cu02 planes. High-quality MgB 2 single crystals were grown using high-pressure cubic anvil technique [33] . 


C. t-J model out of equilibrium 

We investigate the ultrafast energy exchange between doped fermionic carriers (holes) and short-range antiferro¬ 
magnetic fluctuations within the t-J model on a square lattice, 

Htj = -tn Y. + H.c.) -f ^ J(Si • Sj - Jhihj), (1) 

(ij> 

where Ci^£r=Ci,(T(l — is a projected fermion operator, represents nearest neighbour overlap integral, the sum 

(ij) runs over pairs of nearest neighbours and is a projected electron number operator. The 

system is threaded by a time-dependent flux (/)ij(t), which induces the electric held 

Motivated by the experimental results, we model the dynamics far from equilibrium under the following assumptions: 
(i) We mimic the effect of an ultrashort laser pump pulse by instantly increasing the kinetic energy of a doped carrier, 
while the chargeless spin degrees of freedom remain unchanged. In the context of the t-J model in Eq. 0, we 
hrst calculate the ground state with then we perform a sudden change (quench) of the phase at time t=0, 

describing the i5-like pulse of the electric held. Such a phase quench is obtained by setting (t)=7rd(t), where 

ex(y) represents the lattice vector in the a;(j/)-direction. (ii) In a generic situation, promptly after the absorption of 
the pump pulse, the kinetic energy of only a fraction of holes (also referred to as photo-excited holes) is suddenly 
increased. For optimally doped samples (p ^0.16) and the pump pulse with photon energy 2 eV, huence 700 /rJ/cm^ 
and penetration depth 170 nm, the doping rate representing photo-excited holes is rather low, pq ^ 0.01, assuming 
that each photo-excited hole absorbs the entire quantum of light, i.e., 2 eV. We therefore study the dynamics of 
a single photo-excited hole, assuming that the main relaxation mechanism is represented by coupling of the photo- 
excited hole to local antiferromagnetic excitations. As a consequence of finite chemical and photo doping, the average 
number of antiferromagnetic bonds, which are available to absorb the excess kinetic energy of a single photo-excited 
hole, is large but finite. While our method is limited to the calculation of single hole dynamics, we simulate the effect 
of small photo doping by limiting the propagation of the photo-excited hole to a finite L x L plaquette with L < 10. 

Since exact diagonalization (ED) is limited to clusters which are too small to account for the complete relaxation 
within the experimental conditions described above, we extend the lattice size by employing diagonalization in a 
limited functional space. The latter method was recently successfully applied to describe non-equilibrium properties 
of a charge carrier propagating in a local antiferromagnetic background miiii]. The advantage of this method over 
the standard ED in the equilibrium regime follows from a systematic construction of states with distinct configurations 
of local antiferromagnetic excitations in the proximity of the hole. In addition, it remains efficient even when applied 
to non-equilibrium systems on short and intermediate time scales, as long as the antiferromagnetic disturbance caused 
by the local quench remains within the boundaries of generated excitations [mull. We first construct the parent 
state 


= Ck,^|Neel), 


( 2 ) 


which represents a hole with momentum k, doped into the Neel state. The functional space [35] is then generated by 



n-h 


77kin + Hj 




Uh = 0, ...,Nh, 


(3) 


where i7kin represent the kinetic energy term of the t-J model and Hj represents the spin-flip term of the Heisenberg 
part of the t-J model (spin-flip denotes overturned spins with respect to the Neel state; Hj erases two neighboring 
spin-flips created by the propagating hole). The size of the functional space is determined by the parameter Nh- In 
addition, we impose the condition that the largest distance of a local antiferromagnetic excitation away from the photo- 
excited hole, f=(Aj 3 ox, Abox)j should not exceed fVbox = 4, therefore limiting the spatial extent of antiferromagnetic 
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excitations to a 9x9 plaquette. Contrary to the standard ED, we carry out a finite-size scaling not with respect to 
the geometric size of the system but rather with respect to Nh (see Supplementary Information). 

Applying the Lanczos technique we first compute the ground state |'I'(t=0)). The time evolution 
|d>(t))=e“*^*-'*|4'(t=0)) is then implemented using the quenched Hamiltonian. We set the hopping amplitude 
f;i=0.360 eV, which corresponds to the time unit h/th=l-83 fs. At each time step St{th/fi)<^l we generate the 
evolution |\E'(t-i5t))—7'|d>(t)) by using the Lanczos basis [S^. We measure in Figure 2a the time-dependent expec¬ 
tation value of the kinetic energy relative to the initial state, AEkin(t)=Akin(t)-if^°,J, where Ekin(t)=(iLkin(t)) and 
E^°,^=(iJkin(t<0)). In Figure 2b we plot the relative increase of the spin energy on the antiferromagnetic bond 
at a distance r away from the photo-excited hole, AEs(t,r)=Es(t,r)-Es^\r), where ifs(t, r)=J(Si-(.r(t)-Sj+r(t)) 
and £'i°^(r)=J(Si+r(t<0)-Sj+r(t<0)). After the quench, the total energy in the system remains constant, 
AFki„(t)-l-^r AFs(t, r)=const. 


Supplementary Material for: Snapshots of the retarded interaetion of eharge 

earriers with ultrafast fiuetuations in euprates 
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I. SVD ANALYSIS 

The Y-Bi2212uD non-equilibrium reflectivity reported in Fig. lb of the main text, is obtained by combining 

the transient response measured in two spectral windows with laser pulses provided by NOPAl and NOPA2 (see 
Methods). In general, the matrix contains different spectral features and temporal dynamics, which are 

characteristic of the physical processes triggered by the pump beam. In order to analyze the 2D matrix, 

different approaches are feasible. For example, one can model the reflectivity variation at a fixed delay time t, 
i.e., , and describe the entire 2D map as the result of the temporal variation of some parameters used to 

describe the spectral response m- However, this approach turns out to be difficult to apply when the response of 
the sample cannot be described as a single component. In order to determine how many components contribute to 
the full matrix, we adopt a method based on the singular value decomposition (SVD) analysis. This method 

has been previously used to study the interplay between the non-equilibrium dynamics of the pseudogap and the 
superconducting phases in YBi2212 m- It can be shown that a generic m x n matrix, M, can be factorized in the 
form M = TSV where T is a m x m unitary matrix, S is a non negative m x n diagonal matrix which contains 
the singular values (in decreasing order with increasing row number) and V is a n x n matrix. This factorization 
is such that the column of T and the row of V contain, respectively, the trace in the temporal and spectral 
domains corresponding to the singular value of the S matrix. The SVD LAPACK routine is used to compute 
the decomposition and to extract the principal components (singular values) which contribute to the overall optical 
response. It is worth to stress that, since the eigenvalues of the S matrix are non-degenerate, the SVD decomposition 
is unique i.e. it is always possible to find for each singular value a unique couple of temporal and spectral eigenvectors. 
Using SVD, the experimental matrix can be written as: 


^^=^A,vI;^(a;)®<5<(.,(t) (SI) 

eq\ J j 

where kj is the singular value while 'l'j(w) and 5(j>j{t) are the corresponding spectral and temporal eigen¬ 
vectors. The results of the SVD analysis on the Y-Bi2212uD spectral and temporal response are reported in Fig. |S1| 
The main outcome of this analysis is that the first eigenvalue alone accounts for more than 97 % of the global re- 



Probe energy (eV) 


delay (ps) 


0 (rad) 


FIG. SI: a) Transient reflectivity variation of Y-Bi2212uD- b) Time-frequency matrix calculated as the outer 
product between and multiplied by the first singular value fei c) Energy trace at t=40 fs, corre¬ 

sponding to the maximum variation of the transient signal, (orange dots) and the first spectral eigenvector 
(black line) ’I'i(a;) multiplied by (5(^i(t=40 fs) and fci d) Time trace at hoj=l eV displayed with the first temporal 
eigenvector 50i(t) weighted by 4ri(a;=l eV) and ki. For graphical reasons, S(j>i(t) has been offset, e) Weights a, 


of the singular values normalized to one (oy = 
computed for different value of the rotation angle 6. 


E, 


) and sorted in decreasing order f) defined in Eq. 


S2 
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sponse. This result can be appreciated by comparing the raw data with the matrix built as the outer product 
between the first temporal and spectral eigenvectors, weighted by the first singular value = fci'I'i(w) ® S(j}i{t)), 
reported respectively in Fig. |S1^ and b. Both the temporal and spectral behavior is perfectly reproduced by the first 
temporal and energy eigenvector as shown in Fig. |S1| ; and d. The weights of the singular values, normalized to the 
unity, are reported in Fig. Six the weight of the eigenvalues aj, with j > 1, are multiplied by ten to emphasize 


the difference with respect to the first one. In order to highlight the uniqueness of the SVD factorization and the 
effectiveness of the approximation with only the first component, the matrix of the temporal and spectral traces 
are rotated by the same angle 9 and a new matrix ^^(9) is computed as the outer product of the rotated vectors: 
^^{9) = 0 S(f>i{t)R{9) where R{9) is the rotation matrix. The mean-square deviation between the 

experimental data and the computed matrix ^^(9) is: 



(S2) 


and it is reported in Fig. SI' for different values of 9. 


kind of rotation of the eigenvectors, indicating that 


(j{9) displays a minimum for 9=0 and it increases for any 
is the best approximation of the experimental matrix. 


II. TRANSIENT REFLECTIVITY MEASUREMENTS AT LOW TEMPERATURE 

Transient reflectivity measurements on Y-Bi2212uD single crystal are carried out at low temperature in a liquid 
nitrogen cryostat. The minimum temperature achievable in the cryostat is about 100 K. The temperature of the 
sample is monitored by a thermocouple placed on the sample holder. The high temporal resolution (below 20 fs) 
is preserved by using a very thin (200 fim) UV fused silica small flange as optical access to the cryostat. The first 
temporal eigenvectors S(j)i(t) extracted by the SVD analysis of the experimental matrices at different temperatures, 
are reported in Fig. |S2| All the time traces do not display any change in the dynamics in the 100 K-300 K temperature 
range reaching the maximum variation at the same delay time (wdO fs) as shown in the inset. 



delay (ps) 


FIG. S2: Temporal dynamics of eigenvectors at different temperatures. The black lines are the fit to the 

data. The fitting function takes into account the experimental resolution, as described more extensively in 
the main text. In the inset, the dynamics at short delay times clearly show that the rise time of the temporal 
traces does not depend on the temperature. 
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FIG. S 3 : Energy dependence of the scattering rate at T=300 K in the probe spectral range. 700 is the asymptotic 
value of the scattering rate. In the inset the EDM fit to static reflectivity (black line) is compared with 
the calculated reflectivity (red line) under the frequency independent approximation of the damping term: 

liu) = 7 oo. 


III. EXTENDED DRUDE MODEL AND SCATTERING RATE 
A. Equilibrium optical response 

The equilibrium reflectivity i?eq(w), measured by spectroscopic ellipsometry S[3], is reported in Fig. (continuum 
line) of the main text for the Y-Bi2212uD sample (hole doping p=0.13, Tc=83 K). Req{oj) is related to the complex 
equilibrium dielectric function e{u}) by: 


Req (^) 


1 - x/eH 

1 + y/e{u}) 


(S3) 


where = 1 + A model dielectric function, a{uj) = aoioj) + which combines the Extended 

Drude term (see main text) and a sum of Lorentz oscillators centered at energies «1.5, 2, 2.7 eV (crj;,.(w)), that 
account for the interband transitions in the visible range, is used to reproduce the experimental dielectric function 
over the 0-3.5 eV energy range. The bare plasma frequency resulting from the fit is a;p=2.3 eV, in agreement with 
the values from the literature for similar systems close to the optimally doped level S[S]. 

A direct expression of the Extended Drude Model (EDM) scattering rate, 7 ( 0 ;, T), in terms of the electron-boson 
coupling /^y(D) can be derived from the Kubo formula of the optical conductivity in the framework of the electron- 
boson Holstein theory S®. In this case, the electron self energy, S(w,T), can be calculated as a convolution integral 
between J^x(D) and a kernel function L{uj,il,T): 


roc 

E{w,T)= /2x(D)L(w,D,r)dD 
Jo 


where the kernel function can be calculated analytically and it has the following expression: 


L{ijJ,VL]T) = —2Tri 


nin,T) + i 


, 1 D — w 


,1 . D -|- w 

- vp- I - 

2 2ttT 


(S4) 


(S5) 


It is worth to note that the kernel function describes the thermal excitation of both the bosonic modes through the 
Bose-Einstein distribution n(D,T), and the electronic carriers through the terms containing the digamma functions 
dl. Consequently, in the first term of eq. |S5[ T is the temperature of the bosonic excitations while, in the rest of the 
formula, is the temperature of the electronic carriers. Eor moderately doped systems, in which the vertex corrections 
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can be omitted, the optical scattering rate is related to the single-particle self-energy by: 


"f[uj,T) = Im < a; 


^+oo 


m,T)-f{^ + u;,T) 
i + s*(e,r)-s(e + a;,r) 




— CJ 


(S 6 ) 


where / is the Fermi-Dirac distribution, 'S,{uj,T) and E*(w,T) the electron and hole k-space averaged self-energies. 

The bosonic glue, extracted from the fit to the equilibrium dielectric function, extends up to «300 meV, and displays 
a narrow peak at «70 meV. A more detailed description of the spectrum can be found in Refs. S|3l[5l|6]. The 

scattering rate calculated from Eq. ^ is reported in Fig. and displays a variation of about 10 % in the entire probe 
energy window. Because of the weak energy dependency of 7 ( 0 ;), in the 0.7-1.5 eV energy range, the equilibrium 
reflectivity can be reasonably described in this spectral region by the asymptotic value 700 as it is shown in the inset 
of Fig. [S3l 


B. Non-equilibrium optical response and the isosbestic point 

In time-resolved experiments, the pump pulse drives the system out of equilibrium while a second broadband probe 
pulse monitors the transient change of the dielectric function. In general, the relaxation dynamics can follow two 
different relaxation pathways. In the first scenario, the pump energy is absorbed by a small fraction of the charge 
carriers which scatter with the other unexcited carriers and relax on a very short timescale towards a Fermi-Dirac 
distribution characterized by a temperature higher than that the equilibrium. One can rationalize this process in 
term of a sudden increase of the electronic temperature {ST^ > 0 ) while the bosonic bath is essentially unchanged 
{6Th =0). In the second scenario, the excited carriers mainly release the excess energy through scattering with the 
bosonic modes. In terms of the effective temperature models, this picture implies that JTe ~ 5Tb at all the timescales 
longer than the single electron-boson scattering process. The rapid increase of the bosonic excitations results, in turn, 
in the increase of the optical scattering rate (^7oo >0), which can be detected under the form of a broadening of the 
Drude peak. 

As it was already demonstrated (see Fig. I in S[3]), the two physical scenarios previously sketched, result in a 
dramatically different reflectivity variation in the frequency domain. In Fig. I of the manuscript we demonstrate that 
the experimental 5R{uj,i) is reproduced by assuming an increase of 5^co on a extremely short time scale (<50 fs). No 
signature of the electronic carriers decoupled from the bosonic degrees of freedom has been revealed. These results 



FIG. S4: The experimental equilibrium optical conductivity (grey dots) together with the Drude and charge 
transfer peak resulting from the fit of the equilibrium optical conductivity (red and the green lines). The first 
eigenvalue obtained from the SVD analysis is reported (yellow dots). The trace extending between 0.7 and 
1.5 eV is the one previously reported in Fig. 1 of the main text. The trace extending between 1.8 eV and 
2.4 eV results from the SVD analysis of the Y-Bi2212uD transient response probed by NOPA 3. The black 
line is the best fit obtained by increasing the total scattering rate. 
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are compatible only with the second scenario, in which the charge carriers mainly exchange energy with the bosonic 
bath. Therefore, we can assume that 6R{t)oc d7oo(t). 

Having clarified how the transient reflectivity reveals the dynamics of the energy stored in the bosonic bath, we focus 
now on the non-equilibrium optical response around the isosbestic point. In particular, it can be shown that all the 
reflectivity curves corresponding to different scattering rates cross each other at a precise frequency w «1.1 eV, which 
is called isosbestic point. The physical properties of the isosbestic points are discussed in a recent work S[T0]. Across 
an isosbestic point a physical quantity f(u!,g) can be locally approximated by an expansion of a generic parameter 
g around some particular value go- In this specific case, for a small variation of the scattering rate, the following 
expansion of the reflectivity holds: 


OR 

=R{u},-fo) + ■^( 7 - 7 o) + 0 [(7-70)^] 


(S7) 


The effect of the pump pulse, in a pump-probe experiment, is to transiently perturb the optical properties of the 
system. In Fig. 1 of the main text we show that SR{uj,t)/Req{uj) changes sign from negative to positive at the frequency 
corresponding to ui. The major point of the analysis is that the build up dynamics of SR{u!,t)/Req{uj), immediately 
after the pump excitation, can be reproduced solely by assuming a transient increase of the electron-boson scattering 
rate. This result can be quantitatively expressed in terms of the variation of jao- The maximum transient signal, 
occurring at a delay time of ^40 fs for Y-Bi2212uD, is equivalent to a relative variation %. Therefore Eq. 

can be reformulated in this way: 


S7 


6R{uj,t) 1 dR{uj) 
i?e,M ^ Req{uj) 5700 


(S8) 


Comparing Eqs. |S7| and S8 we can easily identify the dynamics of the scattering rate variation S^ao(t) with the 


first temporal eigenvalue obtained through the SVD 64>i{t). 

To further support the interpretation of the experimental results solely as a change of the scattering rate, the 
non-linear optical response is measured in the visible-UV range between 1.8 eV to 2.4 eV using chirped-mirrors 
compressed pulses in a degenerate configuration (both pump and probe pulses are generated by NOPA 3). In this 
energy region the high energy tail of the Drude peak overlaps with the interband optical transitions. In particular, 
a weak additional Lorentz oscillator (LO) at «1.5 eV and a more intense one at ^2 eV are necessary to reproduce 
the equilibrium dielectric function. The peak at 2 eV is usually interpreted as the reminiscence of the charge transfer 
(CT) edge in the undoped compound, arising from the excitation of an hole from the upper Hubbard band to the 
oxygen 2px:^y orbital, and it should be particularly sensitive to rearrangements of the oxygen orbitals. Fig. S4 shows 
that the first eigenfunction (4'i(a;)) of the SVD of the 6R{uj,t)/Req{uj) matrix can be reproduced by assuming an 
ultrafast broadening of the Drude peak, i.e., d7oo>0. This result rules out any possible variation of the position of 
the oxygen orbitals during the ultrafast relaxation process, since the transient variation of the parameters of the CT 
peak would provide a narrow reflectivity variation around 2 eV, incompatible with the measured monotonic variation 
of the reflectivity. 


C. Magnesium diboride 


The MgB 2 transient reflectivity is reported in Fig. |S5^ and is factorized by the SVD. Analogously to Y-Bi2212uD, 
the MgB 2 non-equilibrium optical response is entirely described by the first temporal and spectral eigenvectors as 
shown in Fig. S5 3. The temporal eigenvector S(f>j(t) perfectly overlaps with the experimental temporal traces (Fig. S5 :) 
proving the reliability of the description of the experimental matrix in term of the highest weight eigenvectors. Despite 
its high Tc, it is generally accepted that in MgB 2 , the Cooper pair formation is mediated by optical phonons at 70 meV 
energy, while electronic correlations are expected to play a minor role. The MgB 2 equilibrium reflectivity along the ab 
plane displays a broad plasma edge at a frequency that is blue-shifted («2 eV) as compared to Y-Bi2212uD- Moreover 
an additional Lorentz peak at r:! 2.6 eV is also observed in the visible-UV optical range m due to the cr — 3 tt transition 
close to the M point of the Brillouin zone. Transient modifications of this interband transition are expected to occur 
on an energy scale larger than the probe spectral window and are not expected to influence the non equilibrium 
response. In Fig. |S5^ the MgB 2 equilibrium reflectivity (black line) and the reflectivity obtained by increasing the 
scattering rate (dashed line) in the EDM that reproduces the equilibrium optical data are reported. Also in this 
case, the transient response, reported in Eig. |S5| 3 of the main text, can be rationalized as a broadening of the Drude 
peak, i.e., an increase of the scattering rate. The 7oo(t) curve for MgB 2 shows a build-up time, rr=32±3 fs, that is 
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FIG. S5: a,b) MgB 2 transient reflectivity variation compared with the matrix calculated as the outer product 
between and <501 (t) multiplied by the first singular value k\ c) Time trace at ?ia;=1.3 eV displayed with 

the first temporal eigenvector weighted by ’Fi(a;=1.3 eV) and fci d) Weights aj of the singular values 

normalized to one {aj = ^ ^^2 ) and sorted in decreasing order e) MgB 2 in plane static reflectivity measured at 

T=300 K (black line) and the static reflectivity calculated within the framework of the EDM by simulating 
an increase of the scattering rate (dashed line). 


significantly larger than that observed on cuprates. This physical process is related to the thermalization with the 
optical bond stretching phonons at 70 meV, which display a high coupling strength and are considered at the origin 
of the high Tc of MgBz S[II]. 


IV. TEMPORAL RESOLUTION 


As it was extensively discussed in the manuscript, the build-up dynamics of 7oo(t) is strictly related to the time scale 
on which the charge carriers exchange their excess energy with short range antiferromagnetic excitations. Therefore, 
the experimental temporal resolution is a crucial parameter to obtain a precise measurements of the intrinsic time 
scale of this process. To characterize the temporal resolution of the experiment, the pump and the probe pulses are 
sent through a 20 /rm thick BBO crystal and the spectrum of the sum-frequency generated is detected at different 
delays. The small thickness of the non linear crystal minimizes the temporal broadening and guarantees a broadband 
phase matching for all the spectral components of the two pulses. The cross-correlation frequency resolved optical 
gating (XFROG) trace, i.e. the sum frequency spectrum vs the pump-probe delay, is reported in Fig. S6 The width 
of the cross correlation trace at different wavelengths (Xxfrog) sets the effective temporal resolution of the transient 
response at a given probe wavelength (Xnm hf), according to the energy conservation condition for the sum frequency 


process (. 


-). In Fig. 


S6 we report a single time-trace at fixed wavelength, i.e., Aprobe=930 nm. 
‘P«mp=590 nm, the time-trace at Aprof)e=930 nm has to be compared to 


^ ^XFROG ^probe A 

Since the pump spectrum is centered at A. 
the cross-correlation trace at XxFROG=3b0 nm whose temporal width Atxfrog is estimated to be 19±2 fs. The 
rise-time of the time trace is longer than the convolution (grey line) between the FROG trace (XFROG) and an 
Heaviside function (H), indicating a hnite delay in the maximum variation of SR{uj,t)/Req{uj). 


V. RELAXATION DYNAMICS IN THE t-J MODEL 


A. Finite-size scaling 


We now describe the procedure adopted to extract the relaxation time from real-time numerical simulations of the 
t-J model far from equilibrium. We perform a finite-size scaling with respect to the parameter (see Methods) 
which defines the size of the functional space. For a fixed value of Nh, we extract the relaxation time r(V/,) from the 
time evolution of the kinetic energy A£'kin(t) of a photo-excited hole. We show Ai?kin(t) for different Nh in Fig. S7 1 . 
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Wavelength (nm) 
600 800 1000 



delay (fs) delay (fs) 


FIG. S6: a) Spectrum of the pump (OPA3) and the probe (OPA2) pulses; b) Measured XFROG trace; c) Cutoff 
of the XFROG trace taken at the maximum of the signal (360 nm); d) Transient reflectivity at Aprobe=930 nm 
(blue dots), fit of the data (black line), convolution between Heaviside function H and the XFROG trace at 
I XFROG (grey line). 


Once T{Nh) is obtained, we extrapolate Too = T{Nh —>■ oo) with respect to the scaling with 1/Nh- This procedure 
is analogous to the finite-size scaling in the conventional exact diagonalization where the size of the Hilbert space 
is determined by the lattice size L. In both cases, Hilbert spaces grow exponentially with either Nr or L. In the 
following, we fix the parameter values to the ones discussed in the main part of the manuscript, i.e, we set J/th = 0.3, 
fVbox = 4 (defining the spatial extend of local antiferromagnetic excitations; see Methods) and the initial quench 
energy Ailkin = 2 eV. 

We obtain the relaxation time from two independent and partially complementary procedures. We start with a 
simple exponential fit fi{t) = Aexp (—of the kinetic energy of the photo-excited hole. This is consistent 
with the fitting procedure of the experimental data for Sjao(t), shown in Fig. 1 of the main part of the manuscript. 
However, some care is required when fitting the numerical data with this ansatz. The ansatz is clearly not correct in 
the very first stage of relaxation, i.e., in the time of the order of a single hopping process, t ^ I fs. In addition, in the 
asymptotic limit where the majority of the kinetic energy is already released to local antiferromagnetic excitations. 




delay (fs) delay (fs) 



FIG. S7: a) Kinetic energy of a photo-excited hole, AiJkin(t), for different sizes of the functional space, given 
by the parameter N^- b) The function r(t) = — Ai?kin(t)/Ai5kin(t)) where Ai5kin(t) denotes time-derivative, for 
the same system sizes. Horizontal dashed line indicates the plateau for Nh = 16 where r(t) becomes time- 
independent. c) r vs l/Nh for two different fitting procedures. In Fit 1, we extract T^^\Nh.) from a time- 
independent plateau of r(t), as indicated in panel b). In Fit 2, we extract r^^l(Aj,) from a fitting function 
f 2 {t) = Ao exp {—y/at'^ + +b). Details are described in the main text. Both fitting procedures yield a very 
similar relaxation time, Ttj ~ 15 fs. Dotted lines represent extrapolation for Fit 1, with fitting parameters 
shifted by 3 standard deviations relative to the mean values. 
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a finite functional space limits the spreading of these excitation throu gh t he plaquette, as seen for different values of 
Ai?kin(i) around t ~ 20 fs in Fig. S7 1 . We therefore calculate in Fig. S7 d the function r(t) = — Ai?kin(0/^-^kin(i 
which gives an instantaneous relaxation time. In the case of the exponential function it simply gives r(t) = 

Remarkably, we observe a plateau of T{t) in the time interval 2 fs < t < 10 fs in Fig. |S7[ 3, which indicates that a 
large part of the relaxation process can be to a very good approximation described by a simple exponential function. 
The time-independent regime of t(<) in Fig. S7 d can only be detected for large functional spaces (obtained for 
A/j = 13,14,15,16), which contain over 10® states. Most importantly, the regime of exponential relaxation with 
T{t) = const extends to longer times for larger values of Nh- The values of T^^\Nh) obtained from the plateau of T(t) 
are shown in Fig. S7: (circles). In the latter Figure, the linear extrapolation with respect to 1/Nh provides a very 
accurate fit for the available data and yields = 14.6 fs. 


We complement the htting procedure with a fitting function / 2 (t) = exp (—\/ot^ 
also captures the time evolution of the kinetic energy at the very short time. 


+ b), which, unlike /i(t), 
of the order of t ~ 1 fs. Moreover, it 
obeys the time-reversal symmetry of the investigated quantum-mechanical system S[T^. It gives f 2 it) = 0 at t = 0, 
consistent with our numerical data, as well as it describes a simple exponential decay f 2 {t) exp (—at long 
times, where = Ij^/a. We cut the numerical data at the time t* when T{t) starts to deviate from the plateau 
shown in Fig. [STj o, which depends on Nh and occurs typically between t* G [6 fs, 8 feb The rest of the data is then 
fitted by /2(i)- The resulting values of = l/^/a are shown as diamonds in Fig. S7:. We agai n ex trapolate 
with respect to 1/Nh for large functional spaces which contain over 10® states (solid line in Fig. S7:) and we get 
= 14.1 fs. 


In both fitting procedures described above, the values of and are quantitatively very similar. We therefore 
take Ttj = 15 fs throughout the manuscript to quantify the relaxation time in the t-J model. In addition, even the 
bare relaxation time obtained for a large but finite Hilbert space (r = 9 fs for Nh = 16) yields the characteristic 
relaxation time scale of the order of 10 femtoseconds. This extremely fast relaxation due to the charge-spin interaction 
is the main message of our numerical simulations. 


B. The role of a direct carrier-carrier interaction in the lower Hubbard band 

Motivated by the extremely low density of the photo-excited carriers (shortly carriers) we have so far studied 
a single carrier in the t-J model, thus neglecting carrier-carrier interactions (see the discussion in Methods). In 
this subsection we provide an additional support for this picture and explicitly show that this interaction does not 
influence the ultrafast relaxation of highly excited carriers at least for parameters which are relevant for the present 
experiment. The direct interaction between carriers in the lower Hubbard band is already captured in the original 
t-J model, however it is not obvious how to single out effects originating solely from this interaction. The first- 
guess solution to this problem would be to study how the relaxation dynamics depends on the concentration of 
carriers. Unfortunately such an approach, apart from being technically very demanding, would also be conceptually 
unjustified. Changing the carrier concentration namely modihes the magnetic background that consequently influences 
the coupling between carriers and the magnetic excitations. In order to overcome this problem we have extended the 
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FIG. S8: Kinetic energy Ai5kin(t)/Ai5kin(t=0) in the system of two carriers, for different values of the direct 
carrier-carrier interaction V, see Eq. (S9). The other parameters are the same as in Fig. 2 of the main text, 
i.e., J/th = 0.3, th = 360 meV and Ai?kin(t = 0) = 2 eV. Results have been obtained by exact diagonalization on 
a) 4x4 and b) 6x4 square lattice using periodic boundary conditions. 
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original t-J Hamiltonian by an additional term describing a direct carrier-carrier interaction with a potential V, 


Htj —t Htjv = Htj + V ^ fLifij, 

(hi) 


(S9) 


where the t-J Hamiltonian Htj as well as the electron-number operator fij have been explained in Eq. (1) of Methods. 
After such an extension the interaction between a carrier and magnetic excitations still depends only on J, while the 
carrier-carrier interaction can be tuned by changing either J or V sm- Below we demonstrate that for a fixed J 
the relaxation dynamics barely depends on V, indicating that the strength of the carrier-carrier interaction has no 
significant influence on the relaxation time. It is important to stress that we do not assign any physical relevance to 
the additional interaction term in Eq. (S9). It should only be considered as an auxiliary interaction introduced to 
establish the experimental relevance of results obtained for a single charge carrier in the t-J model. 

Figure [S8] shows the relaxation of the kinetic energy for two charge carriers propagating on 4x4 and 6x4 plaquettes 
solved by exact diagonalization (ED). Note that the carrier concentration in the former case is close to the hole 
concentration in the optimally doped cuprates. We use ED since the two-carrier system on larger plaquettes solved by 
diagonalization in the limited functional space converges much slower than the corresponding single-carrier system. 
As discussed in the Methods, clusters with N < 30 sites solved by ED represent too small systems to account for a 
complete relaxation within the experimental conditions, and therefore the relaxation times should not be compared 
quantitatively to the data shown in Fig. Nevertheless, the main message from the results in Fig. does not 
concern the quantitative results for the relaxation time but it rather tests the sensitivity of the relaxation time on the 
coupling between carriers as modeled by V. We have computed the relaxation dynamics in the presence of attractive 
as well as repulsive potential V with a magnitude that is sufficiently large to overcome the carrier-carrier interaction 
present in the original t-J model. Results are quite independent of V, and different panels in Fig. [^demonstrate that 
the same behavior persists for various sizes of the clusters. We thus conclude that the direct carrier-carrier interaction 
in the lower Hubbard band has a inconsequential influence on the ultrafast primary relaxation in the time-window 
t < 15 fs. However, we expect that this interaction may become important for larger doping (and in particular for 
large fluence), when magnetic correlations diminish while the density of carriers becomes larger. Finally we stress 
that the direct carrier-carrier interaction between particles in the lower Hubbard band should not be confounded with 
the original Hubbard repulsion U in the full (unprojected) Hilbert space, since the latter interaction is actually the 
physical origin behind the coupling between the charge carriers and the antiferromagnetic fluctuations. 


VI. DYNAMICAL MEAN-FIELD THEORY (DMFT) 


The role of short-range AF correlations as the effective bosons mediating the electron dynamics in doped cuprates 
is also supported by the Dynamical Mean Field Theory (DMFT) solution of the Hubbard model S[T4]. Even if an 
accurate characterization of the phase diagram of the copper-oxide planes of the cuprates requires cluster extensions 
of DMFT S[Tn], the single-site DMFT accurately reproduces the evolution of the optical spectral weight as a function 
of both doping s[ini[n] and temperature sm- We consider a two-dimensional Hubbard model with t=0.25 eV, 
t'=-0.05 eV and U=3eV, and we solve the local quantum dynamics using finite-temperature exact diagonalization 
sm approximating the bath with 8 levels. The frequency-dependent optical conductivity is directly obtained from 
the knowledge of the single particle Green’s functions and of the current vertex. Using DMFT, we calculate the 
temperature-dependent optical conductivity, a{uj,T), of a hole-doped material whose properties are assumed to be 
described by a generic single-band Hubbard model. DMFT is used to compute a{uj,T), shown in Fig. S9i, fully 


retaining both the physics of the high-energy transitions at the energy U and that of the short-range antiferromagnetic 
correlations at the energy scale J. At an intermediate hole concentration (p=0.16), cr(a;, T) is characterized by a broad 
Drude-like part {liuj<2A eV) and distinct interband transitions in the 2.4-4.2 eV energy range, that account for the 
localized excitations from the lower to the upper Hubbard band. As shown in Fig.|S9^, the main effect of a temperature 
increase from 250 K (black line) to 400 K (red line) is the overall broadening of the Drude peak with almost no change 
of spectral weight at the energy U. Quantitatively, this is shown by calculating the finite-cutoff spectral weight, i.e.: 


SW{n^,T)= Re<j{uj,T)duj 
Jo 


(SIO) 


as a function of the temperature T. In the simplest version of the Drude model, in which the scattering rate is 
frequency-independent, the approximation: 


swin^T) - ^ 


I - 2 


ijT) 

TtHc 


(Sll) 
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a) b) 
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FIG. S9: a) The optical conductivity of the single band Hubbard Model calculated by DMFT is reported for 
T=250 and 400 K and [7=3 eV. A broad Drude peak extends up to 2.4 eV, whereas the optical proper¬ 
ties are dominated by interband transitions at the U energy scale, b) The relative variation of the optical 
spectral weight, i.e., ASW/SW=SWi(flc,r)-SWi(f2c,250K)/SW(nc,250K), as a function of the temperature T. 
The temperature-dependent intraband spectral weight SWi(2.4 eV,T)=J^''^'^^I{.eai{ij,T)dijj is calculated for the 
DMFT optical conductivity (i=DMFT, red dots) and for an extended Drude model (i=D). The interband 
spectral weight is calculated as the difference between and the intrinsic 

spectral weight of the underlying Drude tail, i.e., (T) =J 2 ^^^J^,jReaDr)cia;. 


is valid in the limit Therefore, the temperature dependance of SW(ftc,T) directly reflects that of j{T). In 

the Fig. SQr, we show that SW(ni=2.4eV,T), calculated from the output of DMFT, is in perfect agreement with 
that expected for a simple Drude model (SW_D(n,T)), provided a suitable effective electron-boson scattering rate 
is i ntrod uced. The quantity SW£)(2.4 eV,T), reported in Fig. [S^ as a red line, has been obtained by replacing in 
Eq. Sll the asymptotic temperature-dependent scattering rate, 7 oo(T) calculated through Eqs. S4 S5 and S6 The 
bosonic function /^x(n) is taken from Ref. [31 The calculated SW_d( 2.4 eV,T) is robust against modifications of the 
details of /^x(n) in the 0-500 meV energy range. Furthermore, the slope of SW(fl2=4.2eV,T) exactly scales with the 
ratio (see Eq. Sll) finally demonstrating that, at sufficiently high temperatures and hole concentrations, the 

dynamics of the charge carriers can be entirely described through an effective electron-boson coupling on the scale 
of J (no other energy scale is present in the model), without any change of spectral weight of high-energy Mott-like 
excitations at U. 
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